On voltage-current characteristics and critical current in Bi-2212 
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V(I) characteristics have been performed in a monocrystalline microbridge of Bi[Ph\ — 2212. The 
vortex phase diagram has been greatly investigated. Linear but non-ohmic Voltage(Current) (V(l)) 
curves with well defined critical current have been observed. A departure from this behavior is 
observed near the peak effect where an out of equilibrium high threshold current can be stabilized. 
At high temperature, the critical current persists in the "liquid" state despite the dissipation at the 
lowest bias. Some implications of these results are discussed. In particular, it is proposed that the 
surface disorder, rather than the bulk disorder, is responsible for the vortex pinning in this sample. 

PACS numbers: 74.25.Sv,74.25.Fy, 74.72. Hs 



The generic vortex phase diagram in cuprates is 
now described in terms of into different thermodynamic 
phases of vortex 0. The main idea is that an or- 
dered Vortex Lattice (VL) is present at low field and 
low temperature, and that it develops into a phase 
possessing a degraded order when the thermal fluctu- 
ations or the static disorder are increased. This cor- 
responds respectively to (ordered solid/liquid) and (or- 
dered solid/disordered solid) transitions. The strong ex- 
perimental facts which have justified these ideas are the 
disappearance of the vortex pinning in the " liquid" phase 
and the increase of this vortex pinning (the peak effect) 
when crossing the disordering transition. The highly 
anisotropic cuprate Bi-2212 can be taken as a representa- 
tive sample where these three states can appear. Never- 
theless, recent experiments suggest to qualify this point 
of view. Indeed, it has been found that both the high 
field and the high temperature properties can be inter- 
preted with the same state of vortex matter, meaning 
that no difference between them should be presupposed 
0, Q • Another striking result is that VL translational 
order is not a good order parameter to characterize the 
first order transition ^ , whereas one could have expected 
the contrary for a genuine melting. From an experimen- 
tal point of view, the VL behavior in Bi-2212 has been 
tackled at low field by magneto-optic imaging or local 
ac probes dQ, which are sensitive to pinning induced 
screening currents on the surface of the samples. If the 
magnetic field is increased, the limited resolution of the 
above-cited experiments strongly restricts their ability to 
give information on the VL pinning and dynamical prop- 
erties. One complementary and easily understandable 
experiment is the measure of a voltage versus current 
(V(I)) curve. The critical current /c, which gives the 
pinning ability of the medium, can be extracted without 
the need of complicated and possibly unjustified assump- 
tions. Furthermore, it has been shown by numerous the- 
ories and simulations 0, IS ^^^^ the nature of the VL 
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should govern the variation of its velocity as function of 
a bulk training force. To the extent that this force is di- 
rectly given by the amount of applied transport current, 
these signatures are directly expected in the V(I) curves. 
For example, there is some consensus concerning the dy- 
namical properties of a VL 0. They are expected to 
be dominated by hopping over barriers at low currents, 
with a transition to a free flow and Ohmic-like regime at 
high current. Fig.l is a schematic representation of the 
different V(I) curves which can be expected for the differ- 
ent VL states, namely ordered state (Bragg Glass) and 
disordered states (vortex Glass and vortex liquid). We 
note that if numerical simulations of velocity versus force 
curves are extremely numerous in the literature, very few 
direct measurements have been performed in cuprates to 
see if the main expectations compare well with the experi- 
ments. In the present experiment, we are thus principally 
interested in the qualitative differences (if any) between 
the functional form of the V(I) curves, for different loca- 
tions (B,T) in the phase diagram. 

The measurement of a V(I) curve is conceptually a very 
simple method and is currently used as a useful probe of 
vortex properties in low Tc materials. On the contrary, it 
was unemployed for high Tc's, except using very low cur- 
rent densities. The reason is practical: quasi-perfect con- 
tacts are usual in low Tc materials {Rcontacts ^ 0.00117) 
but not in cuprates {Rcontacts ^ 1-2 i7). It is thus difficult 
to avoid any overheating in these rather resistive contacts 
as soon as the current exceeds tens of milli Amperes. Typ- 
ically, under normal Helium atmosphere and standard ex- 
perimental conditions, IsslOO mA into a Ifl metallic con- 
tact (< mm^) is compatible with a local increase of the 
temperature of about 0.1 K. This strongly restricts the 
value of the "safe" transport current. Extrapolating to 
a typical sample size (W=l mm, t— 100 fim, W=width, 
t=thickness) , the V(I) curves are thus difficult to perform 
as soon as J > 10^ A/cm'^ {i > Q.5A/cm if expressed in 
surface current units). Since the critical current should 
be at least less than these values, only the depinning on- 
set close to the high temperature first order transition 
can be reliably measured. It is worth noting that sim- 
ilarly low current densities are usually not sufficient to 
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overcome unavoidable critical current inhomogeneities at 
the onset of vortex depinning, even in pure supercon- 
ducting metals. In this regime, only parts of the VL are 
moving. To reach the flux-flow regime where the whole 
VL is in motion, it is necessary to increase the value of 
the injected current, say (/ — Ic) reasonably larger than 
Ic- One of the possibilities which allows to overcome the 
overheating problems is the use of fast current pulses. 
This technique was employed in Both stable linear 
and metastable V(I) curves with a S-shape were mea- 
sured. Since the pulse of current injection is very fast 
(about 100 /is), a normal skin effect can affect the pre- 
ceding results and it is not completely clear if they can be 
taken as representative of the steady state of the moving 
VL. A confirmation of this experiment, using a continu- 
ous transport current, appears thus necessary. For that, 
we have studied monocrystalline bridges of 100 or 200 /xm 
of width. This ensures a good homogeneity of the current 
injection. Above all, for a same current I, a largest cur- 
rent density {I/tW — J or I/2W — i) is obtained when 
using a microbridge rather than a bulk sample. This let 
the possibility of working in the flux-flow regime without 
the overheating at the resistive contact pads, even at low 
temperature. 

In this experiment, we would like to address simple 
experimental facts to the following questions: 

- What are the fundamental differences, if any, in the 
shapes of the V(I) curves for the different VL states? - 
Is it possible to explain the critical current values? 

The samples used in this study are slightly 
Pb doped single crystals of the Bi-2212 family 
{Bii,g,PbQ,2Sr2CaCu20s+5)- They were grown by the 
self- flux technique as previously described 0. Each 
cleaved single crystal was laser tailored in the form of 
a microbridge with a controlled pattern of (W=200 * 
L=400) ^im^ or (W=100 * L=400) ^im^ (Fig.2). The 
thickness is about 100 fim. The crystal was after an- 
nealed under a controlled pure oxygen gas flow and is in 
the slightly overdoped regime {Tc= 79.5 K). Low resis- 
tance electrical contacts (< Ifi) were made by bonding 
gold wires with silver epoxy. The DC transport mea- 
surements were performed using a standard four probe 
method (cryostat Quantum Design with a 9T horizon- 
tal coil, external current source Adret and nanovolmcter 
Keithley). The results presented here correspond to the 
microbridge with W= 200 fim, otherwise it will be spec- 
ified in the text. 

RESULTS IN THE LOW TEMPERATURE REGIME 

We have first performed V(I) curves at low temper- 
ature (T = 5K) in order to minimize the effect of the 
thermal fluctuations. Let us discuss the results for high 
magnetic field values (IT > B < 9T, Fig. 3). The V(I) 
curves are very similar to those observed for a conven- 
tional vortex lattice, as it can be measured in low Tc 
metals or alloys. In particular, they present the usual 
form V = Rf f{I — Ic) as soon as I is slightly higher than 
the critical current Ic {Rff is the flux- flow resistance of 
the VL). There is no evidence of an ohmic regime at low 



applied current and the depinning is rather stiff. Fur- 
thermore, we have tried to compare the effects of a Field 
Cooling (FC) and of a Zero Field Cooling (ZFC), or of a 
FC under different cooling rates. Always the same dissi- 
pation has been measured in the time scale of our exper- 
iment. In particular, no aging effect is observed on the 
critical current. This is not in agreement with a glassy 
nature of the VL governing transport properties. 

When the magnetic field is decreased, a different be- 
havior is observed in a restricted region of the phase di- 
agram (0.05T < B < 0.2T). If the VL is prepared af- 
ter a FC, the V(I) curve exhibits a S-shape with a high 
threshold current Ihigh flO] , but only for the first increase 
of the current. After this initial ramp, a reproducible 
Ic < Ihigh is always measured (Fig. 4). This has been 
previously observed in the pulsed-current experiments, 
and this state with a "high threshold current" has been 
evidenced as a metastable state with a very long relax- 
ation time 1^, ITU . Our measurement using a dc current 
proves that the observation of this state is not due to 
the kind of stimulation used. This observation is also 
in agreement with the observed supercooling of a high 
critical current state using a time-resolved local induc- 
tion measurement 0. We observe a hierarchy in the 
accessible threshold currents. Depending on the exact 
preparation of the FC state (cooling rate, value of the 
initial field, rate of the current injection ramp), numer- 
ous threshold values are accessible between Ihigh and Ic- 
It is clear that in such a regime where it is easy to lock an 
out of equilibrium state, most of the measurements will 
give transient and spurious relaxation effects if the field is 
ramped or quenched as it is done in magnetization mea- 
surements. We find some traces of this metastability up 
to about IT, but the most obvious effects are restricted 
in the region (0.05-0. 2T) (Fig.5). For / > high, the V(I) 
curve is observed disrupted. If the width of the micro- 
bridge is narrower, multiple steps in the current induced 
moving state can even be observed (inset of fig. 4). Even 
small (100 iJ.m), this width is much larger than any super- 
conducting length. It is thus likely that these structures 
in the V(I) curves can not be explained by phase slip- 
page processes as observed in one dimensional supercon- 
ductors ^3 • This mimics rather the voltage steps which 
have been observed in the current-induced resistive state 
of type I superconductor l_3j and can be the counterpart 
of a coexistence between two states possessing different 
critical currents and comparable spatial extensions. This 
dynamic is close to what is observed in a current car- 
rying inhomogeneous superconductor. It can not be ex- 
clude that the high values of the threshold current are 
responsible for local heating in the interfaces between 
the domains and favor the formation of resistive "elec- 
trothermal" structures 0|. Concerning the metastable 
V(I) curves, the peak effect in the critical current, the co- 
existence of two VL states, the same kind of behavior is 
currently observed in 2H — NbSe2 :15]. The strong differ- 
ence is that the peak effect and the associated metastable 
effects appear close to Bc2 in NbSe2 but is here restricted 
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to a very low field value. Since the applied temperature 
is very similar in both experiments, the explanation of 
this field value difference has to be found in the large 
difference between the electronic anisotropics. For a field 
lower than about 0.05T, we do not observe any hystere- 
sis within the V(I) curves. This is summarized in fig. 5 
where Ic versus B is shown. 

One has to remark that the variation of Ic{B), if one 
excepts the small low field part where metastability takes 
place, looks like what is measured in soft low Tc mate- 
rials. To some extent, one can speculate that the same 
pinning mechanism is acting without involving a different 
VL phase. Qualitatively, the functional form of Ic{B) is 
very close to the one of the reversible magnetization of 
a high K anisotropic superconductor |l7l| . meaning that 
Ic is linked to the weight of the diamagnetic screening 
currents. This has to be brought close to the linear V(I) 
curves that we measured. We will return to this point 
later. 

HIGH TEMPERATURE RESULTS: THE "LIQUID" 
STATE 

Let us now discuss the V(I) curves obtained at high 
temperature. The temperature T was fixed at 50 K and 
the magnetic field B was varied from O.OOIT to 9T, in or- 
der to cross the transition between the so-called "solid" 
and "liquid" vortex phases. In Bi-2212, thermodynami- 
cal consistency required by the respect of the Clausius- 
Clapeyron relations (ISj can not be proved, because up 
to now no transition can be detected using specific heat 
measurements. The small step usually observed in the 
magnetization [l^l , added to the reasonable idea that the 
same physics is acting both in Bi-2212 and in YBaCuO, 
can be taken as a good indication that this transition 
is of first order. The appearance of an approximate lin- 
ear resistivity when the sample is probed at very low 
current density is also usually taken as a good criteria. 
The underlying ideas are that a liquid state of 2D vor- 
tices can not be pinned and that the resistive properties 
are close to that of a metal (Ohmic regime). Using the 
same criterion, the transition would be located at about 
0.08 T (fig. 6), in agreement with values currently re- 
ported in the literature (0.065 T in for a Bi-2212 
with an equivalent doping). Nevertheless, we observe 
that for magnetic fields largely higher, in a large part of 
what is identified as the "liquid" state, the V(I) curves 
are clearly non ohmic (fig. 7). A non linear response 
in this resistive state has already been observed. It has 
been interpreted as the feature of a highly viscous liquid 
of vortices lo], or by the edge effect of surface barri- 
ers (SB) j2Qj. Roughly speaking, SB are expected to be 
negligible for fields B > Bcin/ln{K) « 20i?ci for high 
K superconductor (k — 100) with an ideal fiat surface 
[2H. Any real sample also possesses surface irregulari- 
ties which facilitate vortex nucleation and decrease this 
value. We estimate that, taken Bd < 300G, SB effects 
can be neglected when a magnetic field of several Teslas is 
applied. We note also that other experiments performed 
using the Corbino geometry do not evidence any effect of 



SB at high temperature . We have also checked that 
a decrease of the sample width decreases /c, confirming 
that neglecting SB is reasonable, at least as far as crit- 
ical current properties at high field values are involved. 
Wc stress on the following points: the high current part 
of the V(I) curves is linear and its extrapolation never 
goes to zero (fig. 7), at least for B < 9 T. In other 
terms, putting aside the rounded dissipation onset, V(I) 
curves can be expressed asV = Rff{I — Ic) as it was ob- 
served at low temperature. It can not be explained by a 
non linear mechanism in which a depinning energy would 
be a function of the driving force required to overcome 
barriers, because the high current regime is linear but 
does not reach an asymptotic Ohmic regime {V — R.I). 
This demonstrates the existence of a real critical cur- 
rent < /c < 10mA in this high temperature state. This 
critical current exhibits a field dependence that compares 
well with the one obtained at low temperature. It is thus 
reasonable to think that the same pinning mechanism oc- 
curs at low and high temperatures. The difference in the 
low dissipation level can be understandable in terms of an 
additive process which appears at high temperature, and 
a thermally activated process which assists the depinning 
is a natural candidate. This was largely studied in the 
literature 0. This defines a threshold current I* above 
which a small dissipation appears but pinning continues 
to exist (Fig. 6). 

DISCUSSION 

From a theoretical point of view, vortex lattice depin- 
ning is generally described like a critical phenomena: the 
bulk depinning of an elastic system in a random media. 
Driven states are described with an overdamped dynam- 
ical equation and it is expected that without thermal 
activation the velocity v just above the depinning scales 
like [F — Fc)^ where F is the applied force, Fn the bulk 
pinning force and (3 is the critical exponent [l|, [23 ■ The 
usual analysis supposes that [F — Fc)^ identifies with 
{J — Jc)^ , J being a current bulk density. At high ve- 
locity, disorder is found not relevant that leads to a ve- 
locity vaF. Experimentally, this should correspond to 
an Ohmic regime. Clearly, the experimental results are 
different. Even without performing a detailed analysis, 
one can realize that the V(I) curves at the depinning on- 
set are never observed convex. P is never lower than 1 
even at low temperature where the vortex Creep driven 
by thermal activation is negligible. One excludes from 
this discussion the peculiar case of the peak effect where 
a convex part is effectively observed but is very likely due 
to the macroscopic inhomogeneity of the state. It is more 
important to remark that, when I ^ Ic, the velocity is 
found to vary like (I — Ic) and not like /. This means that 
the pinning force felt by the VL remains constant even at 
" high velocity" , in apparent contradiction with an expla- 
nation in terms of a bulk depinning. One could object 
that the applied current in our experiment is not suf- 
ficiently high to reach the predicted Ohmic-like regime. 
This problem can never be strictly resolved because Joule 
heating is always a limitation when increasing the cur- 
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rent in a transport experiment. Nevertheless, one can re- 
fer to low Tj, materials where the experimental situation 
is much more attractive. Under very controlled temper- 
ature, it is possible to verify Va{I — Ic) up to at least 
/ > 30/c. This is thus the generic shape of a V(I) curve 
in a type II superconductor, and the present experiment 
shows that the same dissipation mechanisms are acting 
in Bi-2212 samples. It has been also verified by inelastic 
neutron scatterin g th at, in this regime, the VL is moving 
freely as a whole |2J|. This latter result means that the 
bulk disorder can be estimated to be no more effective. 
At the same time, the former result shows that the ve- 
locity does not reach the asymptotic regime where the 
velocity should vary linearly with the total current. To 
be coherent with the above mentioned theories '23], one 
possible solution is to replace by hand the applied force 
F (resp /) that acts against bulk disorder by {F — Fc) 
(resp [I — Ic)) (Fig. 9). Physically, this solution appears 
in the case of vortex pinning by the surface roughness 
|25| . The main idea is that the surface disorder (a quite 
standard surface roughness) allows, thanks to boundary 
conditions, for the flow of non dissipative superficial cur- 
rent {ic{A/cm) = Ic/IW). This is only when all the 
superficial non dissipative paths are exhausted, precisely 
when I > Ic, that the excess of current (I—Ic) penetrates 
the bulk. To some extent, a bulk force makes then sense 
and a "bulk depinning" can occurs, but involving only 
the over critical part of the applied current. The bulk 
disorder can be averaged by the motion without affect- 
ing the main critical current which reflects the surface 
pinning ability, explaining the experimental shape of the 
V(I) curves. As an further indication that this mecha- 
nism can occur in the Bi-2212 samples, quantitative ex- 
pressions for Ic{B) are predicted and can be checked. The 
case of very anisotropic samples is specially interesting, 
because it is predicted that for not too low magnetic field 
values and for a moderate surface roughness, the surface 
critical current ic may become independent of the surface 
quality and depends then solely on the parameters of the 
condensate |l7l |. For clarity, we restrict the comparison 
to the high field values in order to use the Abrikosov 
limiting expressions p^ : 
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2iiq(3a1k.2//c{T)' 

Here, K2/ /c{T) is the generalized Ginzburg-Landau pa- 
rameter, 7 the electronic anisotropy and Bc2 the second 
critical field. For the small temperature dependence of 
K2, we use the microscopic GL derivation |2y|. This ex- 
pression describes rather well the data (Fig. 10), with the 
use of a restricted number of parameters. 7 = 60 was 
obtained from angular resolved resistivity measurements 
using scaling arguments With K2 = Bc2 = 
we find ^(0) — S.Snm, A(o) = 2080nm consistent with the 
slightly overdoped regime p^). The approximate field 
dependence of the critical current in a restricted field 
range is an indication but can not be taken as a proof of 



the validity of the model. Other field dependence would 
be acceptable. In particular, in a log-log representation, 
it is reasonable to estimate that Ic varies like as 
it is often observed ('s^ and references herein) and inter- 
preted as a 2D strong bulk pinning regime. A possible 
limitation of this latter analysis is the relevance of the 
strong bulk pinning centers in clean single crystals. It 
is more important to note that the absolute value of the 
critical current is the one expected if only the surface 
contributes to pinning. To compare with other experi- 
ments, we measured ic ~ 0.1-0.5 A/cm at (T= 5K, B= 
9 to IT) in our microbridge, and 0.4-1 A/cm was mea- 
sured in 0. These values seem representative of Bi-2212 
crystals. Now it can be objected that the current dis- 
tribution in Bi-2212 can be strongly influenced by the 
very large electronic anisotropy (32j . This leads to a very 
small current penetration in the normal and in the super- 
conducting state and this has a priori no link with any 
pinning mechanism. It remains that if a current flowing 
under the surface is non dissipative up to a critical value, 
a surface pinning mechanism should apply. Finally, it 
seems that there is no need to involve bulk disorder in 
order to explain the critical current and the shape of the 
V(I) curves at least outside the peak effect region. Other 
experiments will soon be performed to evaluate more pre- 
cisely the validity of this hypothesis and the specific case 
of the peak effect will be discussed in another paper. 

It is clear that this interpretation of the VL pinning in 
Bi-2212 is at odds with the currently accepted view that 
the pinning is driven by the bulk disorder (collective bulk 
pinning). We need to add some comments to replace this 
result in the context of recent experiments. 

In Bi-2212, the peak effect should separate a high field 
disordered phase with glassy properties the low field or- 
dered VL. The VL, corresponding to the V(I) curves of 
the fig. 3, is expected to be in this glassy state, whereas we 
observe a conventional VL dynamic rather than a glassy 
dynamic in its proper sense. The experimental proof of 
a disordered high field state in Bi-2212 is also not com- 
pletely obvious. Small Angle Neutron Scattering exper- 
iments can effectively be interpreted in terms of a VL 
disordering transition at very low field p9|. This inter- 
pretation relies in the strong decrease of the diffracted 
intensity when the density of vortices increases, whereas 
a simple London model which disregards the finite size of 
the vortex core predicts only a smooth decrease due to 
geometrical factors. Nevertheless, recent theoretical pro- 
gresses allow to calculate the microscopic field distribu- 
tion and its Fourier components are found to be strongly 
field dependent even at low field (for high kappa and/or 
anisotropic superconductors) ^3Q], meaning that involv- 
ing a transition in the VL may be not necessary. This 
can explain why other experiments with a largest inten- 
sity are showing that a well ordered VL with good Bragg 
peaks survives at higher fields [s^- We note also that 
recent experiments suggest that the VL order is not rele- 
vant for the nature of the "melting" transition in Bi-2212 
meaning that the order parameter of this transition 
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is not linked to the topology of the VL. Experimental 
data suggest that the line of the peak effect in the phase 
diagram is the continuation of the first order line into 
the low-temperature [l^ . This implies that the same 
state of the VL is present at low and high temperature. 
Our experiment suggests that the pinning of the vortex 
lattice is also of the same nature, even if thermal fluctua- 
tions are responsible for a low dissipation background at 
high temperature. Finally, it has been recently proposed 
that the renormalization of the (non-local) line tension 
by thermal fluctuations for large wave vectors is at the 
origin of the first order transition '33| . As the line tension 
is the controlling parameter for the surface pinning, this 
offers a scenario more compatible with our results than 
the occurrence of a genuine melting transition. Further- 
more, even if the bulk condensate is strongly fiuctuating 
with eventually disrupted vortices, a small dissipation- 
less current can persist under the surfaces. A well known 
example is the case of surface superconductivity where 
a critical current does exist in the surface sheath even if 



the bulk develops into its normal state. 

To conclude, the use of monocrystalline microbridges 
of Bi-2212 allowed us to measure V(I) curves in differ- 
ent parts of the VL phase diagram. They appear con- 
ventional and a clear critical current can be defined by 
extrapolating the high current linear part. A metastable 
high threshold state can be stabilized by Field Cooling. 
This corresponds to the peak effect region and mimics 
what is observed much closer to Bc2 in NbSe2. The 
values of the measured critical currents and the shapes 
of the V(I) curves suggest that the pinning by the sur- 
face roughness can not be ignored for understanding the 
vortex lattice dynamic in Bi-2212. We observe that this 
superficial critical current survives in a large part of the 
so-called liquid state. Finally, the overall picture sug- 
gests that the VL pinning and dynamic in Bi-2212 are 
more similar to what is currently observed in low Tc ma- 
terials than it is often suggested. 
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Jl J 2 



FIG. 1: V(I) curves expected for different phases of the VL (A:Quasi-Lattice (Bragg Glass), B:Glass, CiLiquid). Ji represents 
the depinning current in the phase A, J2 in the phase B. 




FIG. 2: The monocrystalline microbridge (200*400 nrri^) of Bi-2212. 
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FIG. 3: Reversible V(I) curves (T = 5 K, B = 1, 3, 5, 9 T from the right to the left). The dashed and straight line is a guide 
for the eyes and evidences the flux — flow regime. 
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FIG. 4: Hystcrctic V(I) curve (T = 5 K, B = 0.15 T after a FC).l/ The first increase of current defines Ih^gh and 2/ The 
following decrease of the current defines Ic ■ In the inset is shown the same kind of curve in the microbridge with a smallest 
width (T = 5 K, B = 0.05 T). Note the steps in the V(I) curve. 
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FIG. 5: Peak effect in the critical current (T=5K. Black points, Ihigh and empty points, 7c). In the inset is shown the variation 
of Ic at high field. 
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FIG. 6: Low dissipation level of the V(I) curves (T=50K, from the left to the right: B= 0.14 T, 0.12 T, 0.1 T, 0.08 T, 0.06 T, 
0.04 T, 0.02 T). As shown in the figure, I* is defined using a criterion of 0.2 fjV and is nearly zero for B > 0.08 T. 




FIG. 7: Reversible V(I) curves (T = 50 K, B = 1, 3, 5, 9 T from the right to the left). The dashed and straight line is a guide 
for the eyes. This linear extrapolation defines Ic- Note that the regime is clearly non ohmic up to at least B = 9 T. 



